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, Abstract. The recent remarkable developments in quantum optics, mesoscopic and cold atom 

(N . physics have given reality to wave functions. It is then interesting to explore the consequences 

of assuming ensembles over the wave functions simply related to the canonical density matrix. In 
this note we analyze a previously introduced distribution over wave functions which naturally arises 
considering the Schrodinger equation as an infinite dimensional dynamical system. In particular, we 
\ discuss the low temperature fluctuations of the quantum expectations of coordinates and momenta 

■ for a particle in a double well potential. Our results may be of interest in the study of chiral 

' molecules. 

o ■ 
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INTRODUCTION 



Quantum statistical mechanics is usually based on the canonical ensemble described 
^ ! by the density matrix = Z^^ exp(— /3H), with Z = trexp(— /3H). In his celebrated 

^ \ book on statistical thermodynamics, Schrodinger Jll] remarked that "this assumption is 

CT' irreconcilable with the very foundations of quantum mechanics" because microsystems 

in general will not be in an energy eigenstate. The reason why this assumption can 
\ be accepted is that a more consistent attitude "leads to the same thermodynamical 

5^ ■ results". The more consistent attitude, according to Schrodinger, consists in introducing 

an ensemble over the states (wave functions) in which one attributes the same probability 
to all the highly degenerate energy states of a macroscopic system. He then shows that 
for a subsystem one expects the canonical ensemble to be approximately correct. In a 
recent paper |0] Schrodinger's statement has been substantially strengthened. See also 
i]. 

In the last decade the remarkable developments in quantum optics, mesoscopic and 
cold atom physics have given reality to wave functions. It is enough to think of the many 
realizations of "Schrodinger cats" or of coherent states. It is then interesting to explore 
the consequences of assuming ensembles over the wave functions simply related to the 
canonical density matrix. A guide in the choice of an ensemble is the observation that 
the Schrodinger equation, idt^f = H\if, and its complex conjugate can be considered as 



^ Dedicated to Giancarlo Ghirardi on the occasion of his 70th birthday. 



an infinite dimensional Hamiltonian system in the variables \]/ and \]/* with Hamiltonian 
(V/,Hi//). It is therefore natural to concentrate our attention on measures invariant under 
this dynamics. In |0] an ensemble over the wave functions was introduced, ^ hereafter 
named Schrodinger-Gibbs (SG) ensemble, with formal measure 

d^SG{¥) = Z^ie-P(^^''^^d{\ - (V/, iir))Yldxir, (1) 

where (•, ■) is the usual scalar product in the Hilbert space. A simple calculation shows 
that in terms of (U]) the canonical ensemble can be written 

Pl3 = J IV^)('^I ^''^^ d^sG{¥), (2) 

where Xj/k are the unnormalized eigenstates of H, assumed to have a discrete spectrum. 

In 10] the motivation was to calculate the distribution of quantum expectation values of 
the operators q and p induced by the SG ensemble. A formula for such a distribution was 
found at low temperature and the wave functions giving the largest contributions were 
characterized in terms of an appropriate Legendre transform of the ground state energy 
of the system in an external field. These wave functions in the case of the harmonic 
oscillator are the usual coherent states and, by analogy, we shall adopt this name also in 
the general case. 

The distribution for the expectation values (q) and (p) is different from that obtained 
with the usual quantum canonical ensemble. For example, in the latter case for reflection 
invariant systems (q) does not fluctuate at all, a rather unphysical result. 

In this paper, af'ter recalling the main steps in 0], we analyze in detail the case of 
a one-dimensional double well, which is an ubiquitous system in physics and deeply 
non classical at low energies. The result is a Gibbsian distribution of the expectation 
values (q) and (p) of the form /'((q), (p)) ^ exp{— /3[(p)^/(2m) + Veff((q))]}, where 
Veff ( (q)) is an effective potential with a single minimum. This minimum is at zero for the 
symmetric double well or near the point where the ground state function is concentrated 
in the asymmetric case. This result may be of interest in connection with pyramidal 
molecules, like ammonia or other molecules potentially chiral. For these systems which 
are adequately described, as far as the inversion degrees of freedom are concerned, by 
a symmetric double well 10,0], fluctuations of (q) correspond to fluctuations of the 
electric dipole moment and could be in principle observable providing thereby a test of 
the statistical assumption (H]). 



LOW TEMPERATURE LIMIT 

We discuss first the harmonic oscillator defined by the Hamiltonian H = p^/2m + 
m(0^q^/2. For this system the Heisenberg equations of motion of the canonical operators 



^ We became aware recently that the same ensemble was introduced also in 101 ■ While some general 
motivations of 0] and jsj are essentially the same, |01 dealt with a specific problem which is further 
analysed in the present paper. 



q and p, namely q = p/m and p = — mco^q, bring to the c-number equations for the 
expectations 



(q) = (p)/m (3) 
(p) = -mw2(q). (4) 

In the following, we will use the notation (q) = q and (p) = p. The Hamiltonian 
flow dHHH) admits the canonical invariant measure Z^^ exp[—l5H{q,p)]dqdp, where 
H{q,p) = p^/2m + m(0^q^/2 and /3 is a constant. We want to show that this result can 
be obtained from the formal SG-measure ([T]). The probability density of the expectation 
values q and p is given by 

P{q.p) = j 5{q-{w.<w))^{p-{w.vw))diisG{w)- (5) 

The Fourier transform of the above expression can be evaluated exactly and transforming 
back to the variables q and p one finds 

P{q,p) = ^exp[-/3(//2m + mft)V/2)]. (6) 

It is interesting to determine which kind of wave functions contribute mostly to P{q, p) 
in the low temperature regime /3 — 0°. To this purpose we have to minimize (v/,Hi//) 
with the three constraints (i//, i//) = 1, (i//, qi//) = q, and (i//,pi//) = p. Once more the 
calculation can be performed exactly and one finds that the minimizing wave functions 
are 



, , /m(0\ 1/4 
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(7) 



i.e. the harmonic oscillator coherent states. 

We wish now to generalize the above results to systems described by Hamiltonians 
H = p^/2m + y(q) with non quadratic potentials V(q). In that case, the density P{q,p) 
cannot be evaluated exactly, however, we can estimate it for low temperatures. In fact, 
for /3 large from steepest descent we have 

P{q,p) - exp[-/3(v/^^p,Hv/^,p)], (8) 

where ^ means asymptotic logarithmic equality and i/z^ p is a normalized state 
which minimizes the expectation of H with the constraints (i//^^, qi//^ p) = q and 
( V9,p; V¥q,p) = P- Fii"st we get rid of the latter constraint by putting 

V^^7,pW =exp(z>;c/r7)(^)^(x) (9) 

with {^q, (l)q) = I. The state 0^ is determined as the ground state ^^^^^ of the eigenvalue 
problem 



{'a + Xq)^X=E^(l)^, 



(10) 



solved self-consistently with the Lagrange multiplier A = X{q) specified by the condi- 
tion 

(<|)J,q<)=^. (11) 
If E^^-j is the eigenvalue associated to the ground state ^i^^y we have 

{v,.p.m,,p) = ^+{^,.m,) = ^+4icj)-H^)^- (12) 

In conclusion, we obtain 

P{q,p) - exp[-/3(pV2m + yeff(^))], (13) 

where the effective potential Veff(^) = ^1(5) ~ ^{^)^ is evaluated by solving the non- 
linear eigenvalue problem (IIOHI 111 . In analogy to the harmonic oscillator case, we call 
¥q.pi^) = 6xp(''P'^/^)<^'A(9) W coherent states. 

On the other hand, in the canonical ensemble for the probability density of the 
expectations of q and p we immediately obtain 

P{q,p)=Z-'l^ e-^^^ 5{q-qk)5{p). (14) 
k 

where Z = Eyt^^^* Qk = {(Pk^^l^Pk)- Here, (pk are the normalized eigenstates of 
H with eigenvalues Ei^. In particular, for system invariant under reflection we have 
P{q,p) = d{q)d{p). 



DOUBLE WELL SYSTEMS 

We consider a particle in a symmetric double well potential of the form 

V{q) = Wo{q^-4)^ (15) 

and solve numerically the nonlinear eigenvalue problem (llOHl II) to obtain the effec- 
tive potential yeff(<?)- This task is accomplished efficiently by the selective relaxation 
algorithm fsl]. The results are reported in Fig. [T] We expect that in this problem one 
can approximate the original Hamiltonian H with a two- state Hamiltonian restricted to 
the lowest two eigenfunctions of H corresponding to the splitting of the ground state 
induced by tunneling. For this reason, we report in Fig.[T]both the exact numerical cal- 
culations and the two-state approximation. In the latter case, the effective potential can 
be evaluated analytically and one finds 



where 



d = (<Phq<P2) 



(17) 
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FIGURE 1. Rescaled effective potential {Ves{q) - (£2 + Ei ) /2) /{{E2 - £1 ) /2) as a function of the 
rescaled expectation q/d evaluated numerically for the double well potential il5\ with different values 
of the mass m. The solid thick line represents the two-state approximation (I16> . For the other parameters 
we set, in all cases, h = I, xq = 1.5, and Wo = 1. The values of Ei, E2, and d used in the rescaling have 
been evaluated numerically and are: Ei = 3.415753, £2 = 4.877688, and d = 1.158335 for m = 0.2, 
El = 2.582908, £2 = 2.865508, and d = 1.268715 for m = 0.5, £1 = 1.970442, £2 = 2.012262, and 
J= 1.353385 form = l,£i = 1.64383345,£2 = 1.65329839, and c/ = 1.38670188 for m = 1.5. 



and ^1 and (f>2 are the lowest eigenstates of H with eigenvalues Ei and E2, respectively. 
Figure [U shows that the two-state approximation is rather good and, in fact, it becomes 
more and more accurate in the semiclassical limit, for example by increasing the value 
of the mass m. 

In the two-state approximation the coherent states i//^ ^ are given by 

xir,,p{x) = e^P^^I'' [ci(^)(pi W +C2(^)(?>2W] , (18) 

where 

Ci{q) = ^-^\-qld C2{q) = ^^/l+q/d. (19) 

In the semiclassical limit the difference E2 — E\ tends to zero exponentially and the 
effective potential becomes flat between ±xq. In this limit J — > xq. This reflects the fact 
that the two levels become equiprobable for low temperatures so that any superposition 
of the corresponding states has the same probability. In the same limit, the dispersion 
— {q^ — q^y/^ tends to xq / V3 as it is apparent in Fig. |2l 
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FIGURE 2. Rescaled fluctuation Aq/d as a function of the rescaled temperature 2/(j3(£'2 — £1)) 
evaluated numerically for the double well potential il5\ with different values of the mass m as in Fig. [2 
The solid thick line is the approximated result obtained from Vgff given by Eq. ilbt . 
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